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Abstract 


We present the numerical simulations of winds from evolved giant stars using a fully non-linear, 
time dependent 2.5-dimensional magnetohydrodynamic (MHD) code. This study extends our 
previous fully non-linear MHD wind simulations to include a broadband frequency spectrum of 
Alfven waves that drive winds from red giant stars. We calculated four Alfven wind models that 
cover the whole range of Alfven wave frequency spectrum to characterize the role of freely 
propagated and reflected Alfven waves in the gravitationally stratified atmosphere of a late-type 
giant star. Our simulations demonstrate that, unlike linear Alfven wave-driven wind models, a 
stellar wind model based on plasma acceleration due to broadband non-linear Alfven waves, can 
consistently reproduce the wide range of observed radial velocity profdes of the winds, their 
terminal velocities and the observed mass loss rates. Comparison of the calculated mass loss rates 
with the empirically determined mass loss rate for a Tau suggests an anisotropic and time- 
dependent nature of stellar winds from evolved giants. 



1. Introduction 


At the latest stages of stellar evolution sun-like stars enter the red giant evolutionary phase losing 
a substantial fraction of their mass in the form of a relatively dense (> 10' 11 M sun /yr) and slow- 
moving (< 100 km/s) wind. Large-scale convective motions of giants and supergiant stars deliver 
CN processed plasma from stellar interiors to the surface, where the wind injects it into 
interstellar medium, where it serves as a source for the formation of new stars and planets. 
However, despite its fundamental importance in the stellar evolution and decades of extensive 
theoretical and observational studies, the physical processes which drive massive and slow winds 
remain unknown. 

Recent high-resolution spectroscopic observations have provided important constraints 
on the physical mechanisms of stellar winds. These observations have revealed four effects which 
are of critical importance in understanding the atmospheric dynamics. First, detailed examination 
of chromospheric emission lines from Fe II, O I and Mg II spectra shows clear evidence for the 
initiation and acceleration of the wind. For giant stars the wind appears to be initiated near the top 
of the chromosphere (Carpenter, Robinson & Judge, 1995). This constraint rules out the 
mechanism of radiative pressure on dust. The mechanism is not capable of lifting atmospheric 
plasma in red giants, because most dust is observed too far away from the stellar photosphere 
(Bester et al. 1996; Danchi et al. 1994). Second, there is evidence for strong turbulence within the 
chromospheres of both giant and supergiant stars. The C II] A.2325 line is formed at temperatures 
between 5,000 and 10,000 K and is optically thin. It has an intrinsically narrow profile, so the line 
width primarily reflects the Doppler broadening caused by atmospheric turbulence. The deduced 
turbulent velocities range from 24 km/s for the K5 giant a Tau to 28 km/s for the K5 hybrid giant 
y Dra to about 30 km/s for the M3. 4 giant y Cru. The fact that the chromospheric sound speed in 
these stars is typically less than 10 km's implies that this turbulence is highly supersonic. This 
constraint suggests that sound waves driven by convection and-'or pulsation do not contribute 



significantly to the initiation of the stellar wind (e. g. Castor 1981; Cuntz 1990; Judge & Cuntz 
1993). While this mechanism can be important in depositing energy in the lower parts of the 
atmospheres of giants and supergiants (for instance, Lobel & Dupree 2001), it also fails to 
reproduce observed terminal velocities and mass loss rates from late-type giants (Sutmann and 
Cuntz 1995). Third, observations suggests that winds from cool evolved stars are at least two 
orders of magnitude cooler than that observed from the Sun and coronal giants (Harper 2001). 
This suggests that the thermal conduction cannot play any important role in driving winds in their 
atmospheres. Fourth, FUSE and HST/STIS observations of red giants have recently revealed that 
“virtually all late-type giants were hybrids to some extent, at least in terms of thelO 5 K 
emissions” (Ayres 2005; Dupree et al. 2005). Observations of hot and broad (>100 km/s) UV 
lines suggests that the effects of magnetic fields can be a crucial factor in dissipating energy in 
the stellar chromosphere and depositing momentum to the stellar wind (Ayres et al. 2003). 

The most promising mechanism to drive massive and slow winds from K and M giant 
and supergiant stars involves incompressible transverse magnetohydrodynamic (MHD) waves (or 
Alfven waves) that propagate along the magnetic field in the stellar atmosphere (Osterbrock 
1961; Hartmann & McGregor 1980; Hartman & Avrett 1984; Charbonneau & MacGregor 1995; 
Airapetian et al. 2000). Dissipation lengths of Alfven waves are much larger compared with the 
compressive acoustic or MHD waves. Early models of Alfven-driven winds from late-type giants 
(Hartman & MacGregor 1 980) suggested that in order to drive massive winds at mass loss rates 


-10" M sun /yr, the total 


enernv flux 




top of the chromosphere should be in the order of ~10 D ergs/cirf/s, where < AFj > is the mean 


square wave amplitude, p 0 and V A0 are the density and the local Alfven speed at the wind 


base. This implies that if the wind base density is p 0 —10' cm" (Robinson & Carpenter 1998), 



then the magnetic fields in the atmospheres of late-type giants should be in the order of -1-10 
Gauss if the area filling factor was included (Harper et al. 2009) . 

This estimate of the magnetic field strength appears also to be consistent with the atmospheric 
magnetic fields required to confine the OV1 emitting plasma observed in cool giants (Ayres et al. 
2003; Dupree et al. 2005; Airapetian & Carpenter 2009). Moreover, GHRS/HST observations of 
non-thermal broadening in Fe II in the spectra of a number of red giants (Robinson et al. 1998; 
Carpenter and Robinson 1999) provide evidence for the existence of substantial supersonic 
turbulence (~ 24-27 km/'s) distributed primarily along the radial and/or tangential directions in 
wind forming regions. The equipartition between the time-averaged (over the wave period ) 


kinetic and magnetic energies of Alfven wave fluctuations, p 0 < SV^ >= 


< sb: > 


magnetic field of ~ 1 Gauss. Recent direct observations of surface magnetic fields in a number of 
giants imply the presence of magnetic fields of a few Gauss (Tarasova 2002; Konstantinova- 
Antova et al. 2009). Another indication of the existence of magnetic fields in cool giants comes 
from VLBI observations of polarized SiO maser emission in evolved stars (Kemball et al. 2009) 
and observations of anisotropic mass losses in from a number of evolved stars as revealed by 
radio interferometric measurements of radial velocities of SiO, H 2 0, and OH masers in AGB stars 
(Hofmann et al. 2000; Thompson et al. 2002). 

Previous models of steady state, spherically symmetric winds driven by Alfven waves suggested 
that in order to produce the low terminal velocities (that are less than half of the surface escape 
velocities) and the high mass loss rates, most of the energy and momentum should be deposited 
by Alfven waves below the sonic point, while momentum addition beyond the sonic point will 
produce the fast and tenuous solar-like winds (Hartmann & McGregor 1980; Holzer & 
McGregor 1985). However, many important questions such as how waves propagate in the 
gravitationally stratified stellar atmosphere and damp their energy were left unspecified because 
they assumed an idealized Wenzel-Kramers-Brillouin (VVKB) approximation. These studies 



parameterized the wave damping and adjusted its level to fit the observations. Later studies, 
including those of An et al. (1990), Velli (1993), Rosner et al. (1995), McGregor & Charbonneau 
(1995) and Charbonneau & McGregor (1995) extended the analysis to non-WKB waves and 
found that waves can be reflected from atmospheric regions where the gradient of the Alfven 
velocity exceeds the wave frequency and that most of the momentum can be deposited below the 
sonic point. These studies used a linear treatment to infer the momentum deposition by Alfven 
waves to the atmospheric plasma to estimate wind acceleration. The location of the sonic point 

(x 

for steady isothermal winds can be expressed as a function of the a parameter as r s = — R, , 

R* 

where R , is the stellar radius, and dimensionless parameter, a = — , characterizes the ratio of the 

H 

kT 

stellar radius to the pressure scale height at the wind base, H = — , g is the gravitational 

R n h,g 

acceleration at the wind base, p is the mean atomic weight of 1.67, and T e is the temperature of a 
gravitationally stratified spherically symmetric hydrostatic atmosphere. This parameter is 
typically greater than 10 for atmospheres of late-type luminous stars. 

To examine the possibility' of producing massive and slow winds by reflected Alfven 
waves below the sonic point, Airapetian et al. (2000) computed the model of stellar winds driven 
by non-linear single-frequency Alfven waves launched at the wind base into the gravitationally 
stratified atmospheres of giants and supergiants using a 2.5D fully non-linear resistive MHD 
code. Our simulations suggested that non-linear effects limit the length of the energy dissipation 
in subsonic regions and provide efficient momentum deposition through the non-linear coupling 
between the azimuthal and radial components of the velocity and magnetic fields in the first 3-5 
stellar radii. This limits the magnitude of the terminal velocity of the simulated wind. To account 
for very large mass loss rates of ~ 1 0 * .V/ / yr that correspond to the wind from a Ori, the 


model required the magnetic field at the wind base, B 0 = 240 G. While such strong surface 



magnetic fields were thought to be unrealistically large, recent MHD dynamo models of late-type 
supergiants suggests the possibility that magnetic fields of comparable magnitude could be 
generated in their atmospheres (Freytag et al. 2002; Dorch 2004). Recently, Suzuki (2007) has 
applied a 1.5D MHD code to simulate winds as well as explain heating rates in the atmospheres 
of evolved giants across the "Linsky-Haisch dividing line" (Linsky & Haisch 1979). His model 
consistently treats acoustic and MHD wave energy dissipation in a fully ionized plasma assuming 
that waves are directly generated at the stellar photosphere by convective motions with the 
magnetic field of ~240 Gauss. However, the photospheric plasma in K-M giants is weakly 
ionized. This suggests due to ion-neutral coupling the damping lengths of Alfven waves are 
about 4 orders of magnitude lower than that expected in fully ionized atmosphere, as shown by 
Khodachenko et al. (2004), and, therefore, Alfven waves will be efficiently dissipated well below 
the top of the stellar chromosphere where the wind acceleration gets initiated. Moreover, because 
the pressure scale height in the cool photospheric layers is at least 4 orders of magnitude smaller 
than the characteristic wavelength of Alfven waves (periods of a day and longer), upward 
propagating Alfven waves will also suffer significant reflection from the regions of Alfven 
velocity gradients. Indeed, if Alfven wave speed jumps discontinuously as waves propagate from 
the stellar photosphere at the Alfven speed, V A .#= 0.7 km/s ( B ph = 1 0 Gs ), to the top of the 

stellar chromosphere at r = 3R, along the radially diverging open magnetic field at the 
local Alfven speed, V M = 138 km's (at B ph = 2Gs), they will suffer the reflection with the 

V -V 

coefficient, R = 1 ’ -0.99 (Davila 1991; Offnan 2002). This estimate suggests that only 

V 4 - V 
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about 1% of the wave flux will leak out to the chromosphere. Therefore, the combined effect of 
the wave dissipation that reduces the wave flux by 4 orders of magnitude and predominant (99%) 
reflection of waves in a weakly ionized and highly stratified plasma of a red giant suggests that 
the Alfven wave energy flux at the top of the chromosphere should be reduced by 6 orders of 



magnitude with respect to the wave flux generated in the photosphere. From the energy flux 
requirements to drive stellar winds discussed earlier, we conclude that the total energy flux at the 


photosphere should be on the order of F — — p n SVlV A ~10 n ergs/cm 2 /s. For the typical 


photospheric density of a red giant and convective velocity ~10 km/s, this requires the Alfven 
velocity, V 4 h ~ 1000 km/s, and, therefore, the existence of the photospheric magnetic fields of 

10 kG! Meantime, the energy flux of photospheric turbulent motions that excite sound 
waves in the stellar atmosphere of a red giant is 



where p h is the photospheric density in a red giant according to MARCS models of 
Gustafsson et al. (2003), V com is the convective velocity. Assuming p h = 3.2 x 10 v 
gr cm -3 according to MARCS models of Gustafsson et al. (2003) and V, onv = 2 km/s 
(Robisnon et al. 1998) in the photopshere a Tau, we obtain the maximum kinetic energy 
flux injected by convective motions ~ 6x10® erg cm' 2 s* 1 . This implies that if the 
compressible waves are converted into incompressible (Alfven) waves in the stellar 
chromosphere with the factor of -0.1%, this can provide the energy flux of Alfven waves 
required to drive stellar winds. Therefore, we assume that Alfven waves that drive stellar 
winds need to be generated at the top of the chromosphere by acoustic shocks produced by sound 
waves formed due to photospheric convection and/or pulsation. 

In this paper we extend our stellar wind models driven by Alfven waves generated at the 
top of the chromosphere to perform a parametric study of the effects of broadband Alfven wave 
frequency spectra on the dynamics of stellar winds, such as terminal velocity and the mass-loss 
rate. Specifically, we develop four Alfven wave-driven wind models that differ by the lowest 



frequency in the Alfven wave spectrum at the wind base and study the trapping efficiency of 
waves and its effect on mass-loss rates and wind speed. We apply our models to K.-M type giants 
stars and discuss their predictions. 

The paper is organized as follows. In Section 2 we present the model equations and boundary 
conditions. In Section 3 we describe the model parameters for the simulated stellar atmosphere. 
The results of our MHD simulations as they apply to the wind from red giants are presented in 
Section 4. Section 5 discusses the implications of the wind models and describes some future 
improvements for wind models. 


2. Time-Dependent 2.5 D MHD Simulations of Winds from a Tau -like Star 

Here we apply a single fluid non-relativistic MHD approximation to treat the dynamics of 
stellar winds extending from the wind base to 20 stellar radii. The validity of this approximation 
is supported by HST/GHRS and STIS observations of dense (n e - 10 9 cm' 2 3 ), slow moving (V < 
100 km/s) stellar wind plasma from K-M giants (Robinson et al. 1998). 



We use a time-dependent single fluid non-relativistic fully non-linear, resistive 2.5D MHD code 
in spherical geometry to simulate stellar wind driven by Alfven waves in the stellar atmosphere. 
We assume that a broadband spectrum of Alfven waves are generated at the lower chromospheric 
boundary and launched along a radially diverging magnetic filed within a specified cone in the 
gravitationally stratified stellar atmosphere of a red giant star. To study the dynamics of wind 
plasma we solve the normalized mass conservation equation (Equation 1), momentum equation 
(Equation 2) and induction equation (Equation 3) for the polytropic equation of state (Equation 4) 
presented as the following 
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dt 
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where the variables have been normalized to the stellar parameters, i.e., radial distance, r = r I R,, 
is expressed in units of the stellar radius, R , , time, t = t It A , in the Alfven transit time over the 

B„ 


stellar radius, r .=R,/V in , where V 


ao 


V 4;r Ai 


is the Alfven speed, and B 0 and p n are the 


magnetic field and the electron density at the wind base, r 1 . The plasma velocity, V = V / V A0 , 
magnetic field. B — B / B 0 , plasma density, p = p! p Q , and plasma pressure, p = p/ p 0 , are 


normalized to V 40 , B 0 , p 0 and the electron plasma pressure, p 0 , at the wind base, respectively. 



The MHD equations (1-4) are written in dimensionless units and contain the following three 
dimensionless parameters: the Froude number is the square of the ratio of the Alfven velocity to 


the escape velocity at the wind base, V 


I2GM. „ . , . . , n . , 

, G is the gravitational constant, R ch is the 

V R ch 


height of the stellar chromosphere, 


K 

vl 


the Euler number is the square of the ratio of 


the sound speed to the Alfven speed at the wind base, E u = 


Pa 

PX o 
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, where c so 


is the 


local sound speed, and the Lundquist number, S = , is the ratio of the resistive 

vc ‘ 

diffusion time to the Alfven transit time and measures the importance of plasma conductivity, 
where v is the classical Spitzer resistivity of plasma parallel to the background magnetic field. 
The maximal value of the Lundquist number is limited by resolution-dpendent numerical 
viscosity in the model. For the grid resolutions used in the simulations presented below we 
simulate the resistive plasma wind at the assumed uniform resistivity corresponding to S= 10 4 . 

The method of solving continuity, momentum and induction equations in fully ionized isothermal 
plasma has been described in detail by Ofman & Davila (1997). The solutions use the 4th order 
Runge-Kutta solver in time with upwind corrections, and 4th order derivatives in space. We have 
also applied fourth order smoothing (i.e., numerical viscosity) to stabilize the solution. The time 
step is determined by the Courant-Friedrichs-Lewy (CFL) condition. 

The equations (Eq. 1) - (Eq. 4) are solved in 2.5 dimensions assuming azimuthal symmetry (i.e., 
5/ 5<p = 0 , two spatial dimensions and three components of velocity and magnetic field). In our 
calculations we assume a uniform, radially diverging background magnetic field to be purely 
radial, the unipolar form 


B o = B o e r / ? ' 2 , 


(5) 



within a cone-shaped atmospheric region with an opening angle of 45 0 (see Figure 1 ) . The radial 
density profile of the initially hydrostatic spherically symmetric isothermal atmosphere is 


Pe = Pn e 


ff(I/r-I) 


( 6 ) 


where a is the normalized gravitational inverse scale height R./H, where HAkT/pm H g is the 
gravitational scale height, g is the gravitational acceleration at the stellar surface, p is the mean 
atomic weight and T is the atmospheric temperature, and k is Boltzmann’s constant. In this study 
we focus on the wind dynamics in red giants driven by Alfven waves on top of Parker’s (1960) 
flow with no specification of heating or cooling mechanisms in the energy equation. We assume 
that the wind plasma is nearly isothermal, i.e., y = 1.01, in the equation of state (Equation 4). 
The boundary conditions at the wind base, r=l, are incoming characteristics approximated by 
zero-order approximation and 

B,.(l,0) = B o , = B„{\,0) = 0, and V r (l,0)>O. 


At the outer boundary, R = 20 R, and at 6 = 0 max we use open boundary conditions with the 
symmetry imposed at 0 = 90° . 

In the present study we generalize our previous wind model (Airapetian et al 2000) by 
introducing the Alfven wave broadband frequency driver at the lower boundary (instead of a 
single frequency driver) as 

B , (t, 0, r = 1) = -V d i V A0 X a, sin (a>,t +r, (0)) (7) 

/'=! 

where a i = i A * , co j = co, + (/ - 1)A®, Aco = (ct) N - co t )/ N - 1, and f; (0) in the phase that is 


the random function in space and time , V d is the Alfven wave initial amplitude. V n) is the local 
Alfven velocity at the wind base, co N = 5co { is the highest frequency in the wave spectrum, and 
/V=100 is the number of modes (Ofman 2004). The wind base in the computational domain is 
divided into 5 regions that generate the Alfven waves with different IT ((?) . 



3. Parameters of Wind Models 


Our MHD simulations of the stellar atmospheric dynamics were performed for an 
axisymmetric single fluid MHD flow in spherical coordinates (r, 0, ip) in the meridional plane 
corresponding to cp^O 0 assuming an initially radially diverging (monopole) magnetic field in the 
gravitationally stratified stellar atmosphere. The calculations were conducted on the two 
dimensional r - 0 grid of 800 uniformly spaced points in radial direction between R=1R« (the wind 
base) and 20/?. and 400 uniformly spaced points in 0 direction covering a cone-shaped 
atmospheric region with an opening angle of 45 0 . 

In this study we simulated the wind from a typical a Tau-like cool giant star, with the model 
parameters presented in Table 1, where the stellar mass, M* is expressed in units of M san , the 
stellar radius, /?. in units of /? sun , the electron plasma density, n e , is in cm' 3 , the sound speed, V s , 
and the Alfven speed, V A0 , are defined at the wind base and expressed in km/s. The models of 
stellar winds were initialized with an isothermal Parker’s (1960) wind flow at T=30,000K. 


Table 1. Summary of the observational and model parameters for the simulated a Tau like 
star 


Parameter 

Initial Value 

Stellar mass, MV M sun 

1.3 (I) 

Stellar radius, RV R su „ 

44 U) 

Mass loss rate ( \-f z i yr ) 

1.6 xlO' 1 "" 

Wind terminal velocity, V, (krn/'s) 

30 (l) 




Mean turbulent velocity, V, (kni/s) 

24 (1) 

Magnetic field at the wind base, B 0 

2 G 

Temperature at the wind base, T 0 

30,000K 

Parameter a 

11.238 

Plasma density at the wind base, n 0 

1 0 9 cm" 1 

Alfven speed at the wind base, V A 

137.99 km/s 

Sound speed at the wind base, V s 

22.25 km/'s 

Alfven transit time, r . 

2.58 days 


Reference: (1) Robinson et al. (1998) 


In the current simulations we use the observational parameters of a typical red giant, a Tau, as 
presented in Table l . This results in a=l 1 .23 at T=30,000K, and the sonic point, r s = 6.62 R, . For 
the given parameters presented in Table 1, the magnetic pressure in the stellar atmosphere 
completely dominates the thermal pressure at the wind base, so plasma motions are directed along 
the radial magnetic field. The radial profile of the Alfven velocity for an isothermal stellar 
atmosphere, with the magnetic field and density profiles given by (5) and (6) respectively, in 
spherical geometry is 


v A 2 e 

r 


( 8 ) 


The left panel of Figure 1 shows that the Alfven speed increases exponentially up to r = — R, or 


about first 3 R , , and then decreases as — 7 - to infinity. The right panel of the figure shows that 

F~ 

the gradient of the Alfven speed reaches the maximum at heights ~2 R, . An et al. (1990) studied 
the propagation of Alfven waves in stratified atmospheres characterized by a parameter as an 
initial value problem and demonstrated that the wave reflection becomes significant at locations 








of large gradient of Alfven velocity. Later, Velli (1993) expanded their model by considering the 
solutions of the stationary problem and studied how the transmission coefficient depends on wave 
frequency. Their analysis suggested that Alfven waves at frequencies, a> ~ (O c , should be trapped 

within the region specified by the maximum of the gradient of Alfven speed. The right panel of 
Figure 1 shows that in the simulated atmosphere of a cool giant Alfven waves should experience 
reflection at frequencies less than the critical frequency, co c = 3.94r^' . For the stellar parameters 

specified in Table 1, the critical frequency is ~17.7|THz with the corresponding wave period of 
4.14 days. 

To investigate the effect of the Alfven wave frequencies on the wind terminal velocity 
and mass loss rates, we constructed 4 wind models driven by broad-band Alfven waves further 
referred to as Models A, B, C and D with the lowest frequencies (and the corresponding wave 
periods) in the simulated broad-band spectrum of Alfven waves presented in Table 2. In Models 
A, B, C and D we assumed the following lowest frequencies of the broad-band Alfven wave 
spectrum, ty, = lr^ 1 , 3 t~ a [ , 9r^' and 18r , or the corresponding wave periods of 16.33d, 
5.44d, 1.8d and 0.9d respectively. This allows us to study the effect of partially reflected waves 
(Models A and B) and freely propagating waves (Models C & D) on the kinematic characteristics 


of stellar winds. 





Figure 1. Initial radial profiles of Alfven velocity (left panel) and the gradient of Alfven velocity 
(right panel) in the simulated atmosphere of a red giant with a=l 1.23. 


Table 2. Wave parameters of wind models 


Model Type 

co | (in r; 1 ) 

Wave period, P (days) 

Model A 

1 

16.2 

Model B 

3 

5.4 

Model C 

9 

1.8 

Model D 

18 

0.9 


The choice of the range of the lowest frequencies, a > s , in the simulated Alfven wave 
spectrum comes from the assumption that Alfven waves are generated by direct perturbations of 
magnetic fields at the wind base driven due to transverse velocity perturbations driven by non- 
linear sound waves propagating upward from the stellar photosphere into the chromosphere. The 
characteristic frequency of the sound wave driver can be estimated from the characteristic 










turnover time of a convective cell. The linear scale of a cell can be estimated from the scaling of 
cells obtained from 3D global hydrodynamic simulations of red giants by Freytag et al. (1997) as 


2.5x10 


-a & M su 

T) nr i / 

*\swi J ej\ sun * 


■R., 


(9) 


where T ejr , and T cfj - stm are the effective temperatures of the star and the sun respectively. 

For convective cell velocities, V ronv ~l0 km/s, the largest convective cell size is l com ~ 0.05 R , , 
which is comparable to the pressure scale height at 7>=30,OOOK. The upper bound for the 

/ 

turbulent turnover time of a convective cell is r = -'-—■■ ■ ■ ~ 1 day . Observations of giants suggest 

' com 

that the frequencies of stellar pulsations in late-type luminous stars is on the order of a few to few 
tens of days. One possible mechanism of Alfven wave excitation at the wind base is due to 
acoustic shock waves that dissipate and heat the stellar chromosphere (Judge and Cunz 1993). 
Another mechanism of Alfven wave excitation that has been suggested is that they are 
parametrically amplified at the half frequency of sound waves, due to a swing wave-wave 
mechanism in the gravitationally stratified stellar atmosphere (Zaqarashvili & Roberts 2002; 
Zaqarashvili & Belvedere 2005). The expected frequencies of generated acoustic waves in red 
giants can also be inferred directly from observations of intensity variations of chromospheric 
emission lines. IUE observations of red giants suggest that the time scale of variations in Mg II 
h&k emission line is ~ 8.5- 12h (Cuntz et al. 1996). Much shorter time scales of emission line 
variations in Ca II h&k emission lines at ~ 1 hour was obtained by Baliunas et al. (1981) from 
power spectrum analysis. Therefore, it is reasonable to expect generation of sound waves in 
giants in a wide range of periods from hours to ~10 days. 


4. Results of MHD Simulations 



We start the simulations by launching Alfven waves with a broad-band spectrum from the wind 
base into the initially hydrostatic spherically symmetric atmosphere with a diverging radial 
magnetic field. The time-dependent solution of Eq. (1-4) is then integrated to a steady-state 
solution. 

Our simulations in Models A through D show that as linear Alfven waves with the amplitude of 
the initial transverse velocity, 8V^= 0.03 V A =4. 14 km's, launched at the wind base propagate 


upward in the gravitationally stratified atmosphere, their amplitude, SV^ , increases with height 
until they become non-linear. This follows from the conservation the Alfven wave energy 


flux,/ 7 = — p 0 SV^V A , for upward propagating non-dissipative waves in gravitationally stratified 


atmosphere, that do not suffer reflection. 

The spatial dependence of the V 9 , V 0 and V r components of velocity and the density in Model D 


in the r-0 plane at t = 192 r A is shown in Figure 2. It is evident that the steady state wind reaches 
the terminal velocity, V m = 0.34r 4 = 46.9 km/s at r~ 10/?, . Figure 3 shows the r-0 profiles 
of the Ohmic dissipation per unit mass, h = j 2 1 p , and the r, 6 and ^components of the local 
normalized Alfven velocity over the computational domain for a steady state wind solution at 
i = 192/^. The left lower panel of Figure 3 shows the formation of localized regions of 
enhanced Ohmic heating at distances greater than 10 /?, from the wind base as the result of a 
non-linear coupling of transverse motions into longitudinal motions in Model D. This is 
consistent with the increase of the Alfven wave amplitudes as they propagate in the 
gravitationally stratified stellar atmosphere reaching the maximum at distances greater than 10/?, 
as shown in Figure 4. The radial profiles of the Ohmic heating integrated over 9 in the simulated 
cone and time-averaged over the longest period of the Alfven w'ave spectrum, P t , 



h = 


Pi 


j -f /(r,g) 


dddt 


in Models C and D are presented in Figure 5. The figure shows that the Ohmic heating,/! , 
slowly increases from distances >5 It in Model D (the model of Alfven-wave spectrum at 
ty, =18 r^ 1 ) as non-linear effects become more important, while in Model C with twice smaller 
<r<>j , wave amplitudes become non-linear only at 15/?,. The heating rate reaches its maximum at 
~15/?»with the values of ,P max =(0.028 -0.04)5^ V J Arc = (1.2-1. 7) xlO 5 ergscnr.T 1 for 

Models D and C respectively. These values imply the formation of warm winds driven by 
outward propagating non-linear Alfven waves. The calculated heating rates in these models are 
about 2 of orders of magnitude greater than that inferred from the surface fluxes observed from 
spectral lines forming in the lower regions of the stellar wind (Carpenter et al. 1999). The actual 
value of the heating rate is proportional to the plasma resistivity, tj . The numerical value of 

resistivity depends on the smallest spatial scale used in our simulations, and is tj - S~ l =10,000 . 
The classical Spitzer resistivity in fully ionized plasma is /7~lCT ll s for the chromospheric 
temperatures of red giants. However, stellar wind plasma in outer regions is partially ionized, and 
the effects of ion-neutral coupling will significantly reduce the plasma conductivity with respect 
to the Spitzer conductivity, and, therefore, correspondingly reduce the value of S by at least 4 
orders of magnitude or to S=10 (Khodachenko et al.2004). In contrast to Model C and D, 
Models A and B with mostly trapped waves output the Ohmic heating that is at least 6 orders of 
magnitudes lower. This suggests that models with fully trapped waves produce much cooler 
winds compared to the models with freely propagating waves. 

We should also note that the energy dissipation observed in our simulations caused by numerical 
compressive viscosity, while the Alfven waves are dissipated by Ohmic heating and phase mixing 



effects because heating is localized also along theta direction can all be important sources of the 
wind heating that will be addressed in the future thermodynamic stellar wind models. 



Figure 2. The spatial dependence of the V ri ,V^ and V r components of the wind velocity normalized to the 
Alfven velocity at the wind base and the plasma density in Model D in the r — 6 plane at t = 1 92? 4 





Figure 3. The spatial r — 6 dependence of the Alfven velocity components, B 0 / -Jp, B t / yfp and 
B r ! tJp , normalized to the Alfven velocity at the wind base, and the Ohmic heating rate per unit mass, 
j 2 / p , normalized to B^ V A / 4/r in Model D at t — 1 92r 4 . 

According to Figure 1, the lowest frequency in the wave spectrum in Model B is close to the 
critical frequency and , therefore, the Alfven waves are expected to be trapped within a few stellar 
radii, while waves at all frequencies in Model D are freely propagating. Figure 6 presents the time 
evolution of the azimuthal velocities in Model B (left panel) and Model D (right panel) as waves 
propagate through the point at r-6.91R, and 0- 90 1 ' (the center of the cone) that 
corresponds to the distance that by 1 R , beyond the sonic point in the stellar atmosphere. 
It is evident that non-linear Alfven waves are observed beyond the sonic point as the 
wind reaches the steady state in Model D, while most of the Alfven wave flux in Model B 




is trapped below the sonic point well before the wind becomes steady state. The 
calculation of the total wave flux in Models A and B shows that only less than 0.1% of the total 
energy penetrates distances greater than 1 R, from the wind base. 



Figure 4. The radial profile of V, component of the steady state wind velocity in Model D at t = 1 92 t A for 
the center of the simulated cone (solid line) and 0 = 45° (dashed line). 




Figure 5. The radial profiles of the Ohmic heating pe 
steady state wind in Model C (solid line) and in Model J 




Time in Alfv#n times Tim® ir. Aifvsn times 

Figure 6. The time evolution of the azimuthal component of the wind velocity in Model B (left panel, 
reflected waves) and Model D (right panel, propagating waves) at r=6.97R* , 0 = 90°. 


Figure 7 shows the radial profiles of the wind terminal velocities in Models A, B, C and D. The 
figure shows that radial acceleration of stellar winds in the first 5 stellar radii is generally 
insensitive to the lowest frequency of Alfven waves, while the wind terminal velocity increases 
by 20% in Model D with respect to Model A. The figure shows that in the models with trapped 
waves (Model A and B), the wind radial velocity is about 30 km/s at 10/?*, while the wind 

terminal velocity -42-45 km/s at 20 R, . The Models C and D outputs greater wind velocities 
reaching over 50 km/s at 20 1L. The calculated velocities at 107?, in Models A and B are 
consistent with the terminal wind velocity for a Tau like star (see Table 1, Robinson et al. 1998). 




Figure 7. The radial profiles of radial velocities of the steady state winds in Model A (solid), Model B 
(dash-dot), Model C (dot) and Model D (dash-double dot). 




Figure 8. The total mass loss rate as a function of the lowest frequency in the broad-band Alfven wave 
spectrum (Model A, B, C and D). The figure shows the sharp increase in the mass loss rate in the model 
with mostly trapped waves. 


The time-averaged, normalized mass flux per unit surface area of the steady state wind is 

WV, =- 7 - }p(t)V r (t)dt, (10) 

where integration of the mass flux is performed over 10 wave periods. The Figure 8 shows the 
total mass flux, M 0 , as a function of the lowest frequency in four simulated broad-band spectra 





(Models A through D). The figure shows that as waves get reflected at lower regions in the stellar 
atmosphere where density is greater, they drive more massive winds. Summarizing our findings 
from Figures 7 and 8, we conclude that our simulations output faster and more tenuous (lower 
mass-loss rate) winds as the lowest frequency in the spectrum of waves increases and waves 
become less trapped in the stellar atmosphere. Our results derived from fully non-linear 
simulations, therefore, are consistent with earlier studies (e.g. Holzer & McGregor 1985) . 

The normalized mass flux per unit surface area, M 0 , in Model A (the wind model driven by 

fully trapped Alfven waves) is a factor of ~3 greater than that calculated in Model D (the model 
of freely propagating waves). The total mass loss rate from the stellar surface in the star filled by 
an open magnetic field is the given as 

M =M 0 

where, <D 0 = 4 n jB^R: , is the total magnetic flux of the region(s) in the star filled with 
the open magnetic field with the filling factor / and the magnetic field, B 0 , at the wind 
base, p 0 is the density at the wind base. Comparison of the calculated mass loss rate from 
Equation (10) with the one observed in a Tau, M ohs = 1.6x 10 u yr (see Table 1), 
yields the average value of the open magnetic field over the stellar surface, the fB 0 ~ 0.75-0.8 

depending on the trapping efficiency of the Alfven waves. For the assumed value of the radial 
magnetic field of 2Gs, f~ 0.37-0.4, that implies an anisotropic wind from a cool giant star like a 
Tau. Thus, we conclude that Model A and Model B output the wind velocities and the mass loss 
rates consistent with the observed parameters for an a Tau-like red giant. 

Moreover, Equation (11) also suggests that the mass loss rate may vary on timescales of the 
emergence and evolution of the magnetic flux associated with an open magnetic field or on 
timescale of the rotation period of the star that reveals changes in distribution of an open 




magnetic flux of the active regions in the stellar atmosphere. While the timescale of magnetic 
field amplification by giant cell convection is in the order of 25 years (Dorch 2004), the rotation 
periods of a typical giant (like a Tau) is on the order of 1 year. This suggests that asymmetries in 
the filling factors of open magnetic filed in two hemispheres can explain variations of the 
measured mass loss rate as ~ <t> 0 , as the star rotates about its axis. A possible indication on the 

role of variations in the filling with stellar rotation comes from HST/GHRS observations of Fe II 
and CI1] lines forming in the wind of a Tau that have been taken in 1990, 1994 and 1995 by 
Robinson et al. (1998). For instance, while the line profile CII] line at three epochs did not 
change, the integrated line fluxes in 1994 and 1996 increased by 1 0%-30% with respect to 1990 
with variation of the UV continuum (2320-2360A) emission in 1994 to 1996 by a factor of 2 and 
1 .77 respectively. Plasma density derived from the optically thin and density' sensitive ratios of 
the CII] lines at three epochs shows no change at three epochs. If we assume that variations in 
line fluxes observed at three epochs were caused by variation in heating rates due to Alfven wave 
dissipation, then the total line flux over the stellar surface, A, is 

F ,=\ jPo < SV; > V ,M rx P,<5V;> <D 0 (12) 

L A 

Equation (12) then suggests that increase in the line fluxes is due to the increased magnetic flux at 
three observed epochs by 10%-30% in 1994 and 1996 relative to 1990 because the density and 
the average amplitude of unresolved wave motions that contributes to the non-thermal broadening 
remain the same. 

The variation of magnetic flux with stellar rotation and/or stellar magnetic cycle that follows from 
Eq. (11) also can explain mass loss rate variations in X Vel observed at 3 epochs, a factor of 6 
variation in the mass loss rates over 18 months in K757 (Ml 5) and by a factor of 2 in K.341 
(M15) and L72 (M13) (Mullan et al. 1998; Meszaros et al. 2009). 



4. Conclusions 


In the current study we simulated slow and massive winds from evolved giants driven by 
non-linear broadband Alfven waves propagating in the gravitationally stratified stellar 
atmosphere of a red giant with stellar parameters of a Tau. The simulations were performed with 
our 2.5D time-dependent fully nonlinear MHD code for four time-dependent wind models 
characterized by the lowest frequency in the broadband spectrum of Alfven-wave driven winds. 
We found that fully reflected non-linear Alfven waves deposit energy and momentum well below 
the sonic point ( Model A & B) and produce ~3 times more massive and 20% slower wind than 
that formed by the fully propagating Alfven waves (Model C & D). This result is consistent with 
solar wind studies and our previous studies of stellar winds (Ofman and Davila 1998; Airapetian 
et al. 2000) and, for the first time, provides the estimate of the contribution of trapped Alfven 
waves to the mass-loss rates from red giants. The calculated mass-loss rates and terminal 
velocities consistently reproduce the observed wind velocity and mass-loss rate from cool giants 
similar to a Tau. Specifically, when the simulated total mass loss rate is compared with the 
measured mass loss rate in a Tau, we conclude that the fB factor, the average value of the open 
magnetic field over the stellar surface, points at an anisotropic wind from the star with the filling 
factor between 0.37-04. In addition, our simulations show that the mass loss rates from Alfven 
wave driven winds depend on the total magnetic flux in open magnetic field in the stellar 
atmosphere that can be either rotationally modulated or produced by dynamo motion and 
transport to the stellar surface. Therefore, our model predicts that the stellar winds driven by non- 
linear Alfven waves from red giants are expected to be anisotropic in space and variable in time 
on a timescale of stellar rotation and stellar cycle. The results obtained in this paper are also 
consistent with detections of anisotropic mass-loss from M giants using diffraction-limited IR 
imaging at Keck (Tuthill et al. 1998: Hofmann et al. 2000) and observations of time varying mass 
loss rates from a number of late-type giants discussed in the paper. 



The current simulations represent dynamic properties of 2D nearly isothermal winds. However, in 
the near future, we are planning to include the effects of Alfven wave dissipation in the energy 
equation due to combination of Ohmic, and viscous heating mechanisms and wave phase mixing. 
Thus, the addition of the energy equation with improved equation of state will be capable of 
reproducing not only dynamics but also thermodynamic parameters of winds from red giant stars. 
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